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We show using diagramtic arguments that in some (but not ah) cases, the temperature 
dependent part of the chiral vortical effect coefficient is independent of the coupling constant. 
An interpretation of this result in terms of quantization in the effective 3 dimensional Chern- 
Simons theory is also given. In the language of 3D, dimensionally reduced theory, the value of 
the chiral vortical coefhcient is related to the formula = —1/12. We also show that in the 

presence of dynamical gauge fields, the CVE coefficient is not protected from renormalization, 
even in the large N limit. 

1 Introduction 

Recently, the manifestation of chiral anomaly in hydrodynamics has been under considerable at- 
tention. Since the original observations in gauge/gravity duality [1-3], it has been shown that the 
hydrodynamic equations of systems with anomalies necessarily contain additional terms if both 
the equations of anomalies and the second law of thermodynamics are to be preserved [4]. These 
new terms lead to two new effects: the chiral magnetic effect (CME) and the chiral vortical effect 
(CVE), i.e. the appearance of a current in the presence of a magnetic field or a vorticity of a fluid 
flow. These effects are dissipationless and hence are proportional to equilibrium thermodynamic 
quantities rather than kinetic coefficients [5,6]. 

It has been found that anomaly matching within hydrodynamics fixes the coefficents of the 
CME and CVE up to a single numerical coefficient [4,7]. This coefficient determines the magnitude 
of the CVE at zero chemical potential. Consider a fluid flow with a velocity profile u^{x). The 
coefficient under consideration is Ci in the following relationship 



Landsteiner et al. [8, 9] show that in theories with gravity dual, the value of this coefficient in 
holography (i.e. at strong coupling) coincides with its value at zero coupling, suggesting coupling- 
constant independence (For a different approach based on kinetic theory see [10]). In a theory 
with a single fermion of right-handed chirality, Ci = 1/12. Landsteiner et al. suggested that this 
coefficient is related to a gravitational anomaly (more precisely, the gravitational contribution to 
the divergence of an axial current). Although the connection between the zero-chemical-potential 
CVE and the gravitational anomaly is direct in holography, it has not been proven outside of 
holography. In fact, naive power counting indicates that the effects of gravitational anomaly can 
appear only in higher-order (namely, third-order) hydrodynamics. 

In this paper, we demonstrate that in a system of fermions coupled through Yukawa interactions, 
the CVE coefficient at zero chemical potential is independent of the coupling-constant. We first 
provide a perturbative analysis of diagrams to show that the one-loop fermion graph is the only 
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graph contributing to the relevant Green function appearing in the Kubo's formula for zero-// 
CVE coefficient. We then show a connection between the CVE and the 3D Chern-Simons term 
appearing in the dimensionally reduced effective field theory of the thermal system. Using this 
connection, we show that we can understand the numerical value of C from the summation formula 
n = —1/12. (The non-rcnormalization in this case can also be inferred from the results 
in [11], where a different approach based on fiuid dynamics in terms of group- valued variables was 
taken.) 

However, it turns out that for the case of coupling to gauge fields, the CVE coefficient is not in 
fact protected against renormalization. Wc demonstrate that even in the large limit, there is a 
single class of diagrams that contribute to the CVE coefficienct and explicitly evaluate this at the 
two loop order. 

The paper is organized as follows. In Sec. 2 we present a perturbative proof of the non- 
renornialization of the chiral vortical coefficient Ci. Our proof is based on similar proof, by Coleman 
and Hill, of the non-normalization of the fermion-induced Chern-Simons terms in (2-|-l) dimensions. 
We elaborate more on the connection to the connection to the 3D Chern-Simons theory in Sec. 3. 
We conclude in Sec. 4. 

Note: As this work was nearing completion, we were informed of [12] which deals with issues 
similar to the ones considered in this paper. The first version of this paper included an error 
which was kindly pointed out by H. Ren et al. who also have a forthcoming paper discussing its 
implications [13]. 



2 Calculation of the CVE coefficient 

In this section, we look at two cases. First we consider a Dirac fermion coupled to scalars through 
a chiral Yukawa term and show that the CVE is independent of this ineraction. In the second part 
of his section, we show that this is also the case for fermions coupled to a dynamical gauge field in 
leading order in the large N expansion, where N is the number of colors associated with the gauge 
field. 



2.1 Coupling To Scaleirs 

To be concrete, we first concentrate on a simplest interacting quantum field theory with a conserved 
axial current; namely, a linear sigma model which contains a Dirac fermion interacting with a scalar 
field. 
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Wc have include in the action the couphng to external axial U(l) gauge field and small metric 

perturbation /i^i^. 

The symmetries associated with the background fields are anomalous, and the anomalies are 
captured by the Ward identities: 



9,4 = - 7^^"'"'F,'^Fp'r, (4) 
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g^j.,. =ir-P( _ ^^r^A^F^i,) , (5) 

where Jg* = V>7^7^V^, F/^^^ is the field strength associated with We note that the axial current 
is conserved in the absence of sources. 

The goal here is to calculate the chiral vortical coefficient at zero chemical potential (^o = 0). 
As explained in Refs. [9], this coefficient is not a proper kinetic coefficient, but is basically an 
equilibrium quantity given by the behavior of the retarded two-point Green's function between the 
current J\ and the momentum density T^^ at zero frequency and small momenta, 

G"^' {u;,k) = ieijnKal + 0{k'') (6) 

ui=0 

The fact that we are interested in the zero-frequency limit allows the chiral vortical coefficient to 
be related to an Euclidean Green's function 

iG^^ = g''^^ = — — Z= {J\T^^) + contact terms, (7) 

oAi ogoj 

where Q^'^^ is the Euclidean (Matsubara) Green's function. At zero frequency we can take our 
source fields and h/^i, to be time independent from the start. With this configuration, it is easy 
to see that the right hand side of the anomalous Ward identities in equations (5) vanish and the 
symmetry is effectively restored. 

We will demonstrate diagramatically, using an argument similar to that of Coleman and Hill 
who showed the non-normalization of the Chern-Simons term obtained by integrating out massive 
fermions in 3D^, that only diagrams that are present in the zero-coupling limit give non-zero 
contributions to the CVE coefficient [15]. The calculation at zero coupling is given in [8], and our 
argument would show that the calculation is exact at any value of the coupling constant. 

Let us now look at a generic diagram at non-zero coupling (see figure la). All such diagrams 
can be obtained from the n-scalar effective vertices with exactly one insertion of Jg* and T^i, by 
contracting external scalar potentials and integrating over the scalar momenta. Let us look at the 
central block of the construction: 

^if{P,Q,Kl,--- ,l^n)- (8) 

Graphically, this is a one-loop graph with n external dynamical scalars, 1 external graviton 6j and 
1 external source gauge boson ai (figure lb). 



^It is known that there are caveats to the Coleman-Hill theorem (see for example [14]). Namely when there are 
parity-violating interactions in the initial Lagrangian. However, as there are no such terms, they do not concern the 
treatment in this paper. 



3 




(a) A generic diagram. (b) n-scalar effective vertex. 



Figure 1: (a) A generic diagram at finite coupling. Tlie internal lines are the sccalar field and the 
external lines are the graviton and the axial U{1) field, (b) An n scalar effective vertex. We can 
get all diagrams of type (a) by contracting the external scalar lines on type (b) diagrams. 



The gauge invar iance of the time independent source fields implies: 

P''^tj\p,Q,l^l,--- ,Kn) =0, 

Q^'^tj\p,Q,i^i,--- ,Kn) =0. (9) 

Differentiating these equations with respect to p^' and and then letting p and q be zero respectively 
we get: 

'^") = ^ rg.i...fe„(p,^,'^i,--- ,Kn) = ^^^^ 

where we have used the fact that these functions are analytic in k. For n = the two momenta p 
and q are not independent and hence we derive: 

r\%,...,Jp,-p) = Oip). (11) 

However, for n > these two momenta are independent and (10) implies: 

rS,..fc„(p,g,'^i,--- ,k„) = o(m). (12) 

This is all we need for our proof. The full contribution to the CVE coefficient is the sum of the 
zero coupling diagram calculated above and diagrams with internal scalars. The latter can be 
constructed from our effective n point vertices by contracting the external scalar lines. 

Two cases arise. If the source fields Oj and 5j are attached to the same fermion loop, the total 
diagram is 0{k'^) by equation (12). Using the same argument, if they are attached to two different 
fermion loops, we again get that the total diagram is 0{k'^) by applying equation (10) for each loop. 
The same reasoning can be applied to the case with the scalar stress tensor insertion and the case 
with scalar internal loops, implying that all such diagrams are 0{k'^). On the other hand the chiral 
vortical coefficient is the term linear in the external momentum of the two-point Green's function. 
Thus, only the one-loop free diagram can contribute to the CVE. This completes the proof. 
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It is important to note that the analyticity arguments given above require the mass of the 
scalars to be non-zero. Given the fact that scalars genericaUy acquire mass at finite temperature, 
this requirement is equivalent to the statment that we are not sitting at a second order phase 
transition. 

2.2 Coupling To Gauge Fields 

We now turn to the case with fermions coupled to dynamical gauge fields. As claimed in the 
introduction, in this case, the CVE in fact does receive corrections. We show that for the non-abelian 
theory at leading order in there is a single class of diagrams that give non-zero contributions 
to the CVE coefficient. To demonstrate this, we point out where the arguments of the previous 
section are affected. The action is: 

S = Jdt d^x {i^0i; - 4^ V^'^" + eA^J^ - ^ V^"*" + > (13) 

where D and V^i, are respectively the covariant (vector-like) coupling and curvature tensor associ- 
ated with the dynamical gauge field V^. We assume one flavor of fermions and we suppress the Lie 
group indices. 

Following the steps of the scalar-coupling case discussed above, we come across two issues and we 
tackle them one by one, one of which ultimately leads to radiative corrections. The first difference 
is that there are now potentially gaplcss excitations associated with the gluons. The infrared 
singularities associated with such massless modes threaten to ruin the analyticity arguments that 
were crucial in the proof. This is the same caveat that also came up in the original Coleman and 
Hill paper [15]. However, since we are working with non-abelian theories, the solution provided 
there does not apply. 

Fortunately, non-abelian theories at finite temperature do produce a mass gap non-perturbatively 
[16]. To take advantage of this fact, we divide the arguments of the previous section into two steps. 
First, we perturbativcly integrate out the fcrmionic fields to give effective vertices with gauge fields, 
the axial source and the graviton (respectively V^, Ai and fiQi) as external lines. According to the 
previous section, the diagrams with more than two external legs are 0{k^). 

Next, we evaluate the expectation value of these effective vertices non-perturbatively. This 
method allows us to use both our pcrturbative diagramatic arguments of the previous section while 
steering clear of the IR singularities that plague perturbative analysis of gauge fields. 

The second difference is that the anomaly (equation (5)) no longer vanishes in the absence of 
external sources and is now given by: 

9,4 = -iCf^^^^'^^M^V + ■■■, (14) 

where we have ommited the terms higher order in the external sources. This is a big problem as 
the arguments of the previous section all depended on the conservation of the axial current. 

At this point, we split the diagrams into two groups depending on the form of the Energy- 
Momentum tensor (see figure 2). The first group has a fermionic Energy Momentum insertion 
(figure 2a) and the second has the gauge part (figure 2b). For the first group, at leading order in 
1/N, we have a single fermion loop to which both the axial current Jg and the energy-momentum 
tensor T"* arc attached. The anomalous contribution of such one loop diagrams is exactly captured 
in equation (14). Expanding the right hand side in metric perturbations, we see that to linear 
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(a) Fermionic stress tensor diagrams. 



(b) Gauge stress tensor diagrams. 



Figure 2: The two diagram classes. 



order, the anomaly only depends on the trace of the perturbation and thus does not contribute to 
correlators involving a single insertion of T*^* . 

However, it turns out that there are certain diagrams in the second class that do in fact con- 
tribute to the CVE. Since in the absence of the anomaly in equation (14), there would be no 
radiative corrections, it is clear that the only contributing diagram are those that include a triangle 
sub diagram as in figure 3a. Figure 3b shows the leading order diagram in this class which we now 
proceed to calculate. 

Since we already know that the only contribution comes from the anomaly in (14), we can 
replace the triangle part of the two loop diagram with the effective vertex: 



Anom 



47r2, 



-9 



(15) 



which captures the divergence of the axial current, and we are effectively left with only a one loop 
diagram (figure 3c) . The calculation of the diagram is slightly complicated by the fact that different 
components of the gauge field receive different effective masses which are in general very important 
for our arguments. However, in this case it turns out that there are no infrared divergences even 
at zero effective mass and we present the calculation with this simplification. 




(a) Anomalous diagrams. 



(b) Leading order diagram. 



(c) Effective leading diagram. 



Figure 3: (a) The only class of diagrams with non- vanishing corrections to the CVE in the large 
N limit, (b) The leading order diagram, (c) After replacing the leading order diagram with the 
effective vertex JAnom- '^^^ triangle vertex represents the anomalous effect of the triangle diagram. 



The anomalous contribution to the Euclidean Green's function is: 



EL 

3 



ni,oj _ / Ti rpOjK _ -9 TC{r)d{G) 



+ 



(16) 



where d{G) is the dimension of the adjoint representation and we have expanded up to first order 
in k. It is now a simple matter to carry out the sum and the integral: 



Anom 



-9'C{r)d{G) / , 1-e- 



2-2- J ^ V 

- Act E 



9^C{r)d{G) 
487r4 



Jjmi. 



(17) 
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where in the last step we have removed the divergent zero temperature contribution before doing 
the integral. In fact, at zero temperature, this contribution is not allowed by full Lorentz symmetry 
and vanishes (This is also verified by direct calculation) . Therefore, any regulator that respects full 
Lorentz symmetry in the zero temperature limit would reproduce the same result. 

We now read off the leading order correction to the CVE coefficient in the large N limit: 

V _ g^C{r)d{G) ^ 

"Correction — ' ^ ' 

3 Connection To 3 Dimensional Chern-Simons Term 

The diagramatic arguments put forward in section 2 were originally presented as proof of non- 
renormalization of the U{1) Chern-Simons term in 3 dimensions. It is well known that the topolog- 
ical nature of the non-abelian Chern-Simons term leads to its quantization and therefore, does not 
receive radiative corrections. Similar arguments for the Abelian case can be made if one assumes a 
compact spatial manifold. 

In this section, we will show that under similar assumptions, one can show that the quantization 
of the Abelian Chern-Simons term in 3-dimensions leads to the non-renormalization of the CVE 
coefficient. First, we show that gauge invariancc leads to the quantization of the level and then, as 
an example, we verify that the free theory satisfies these quantization conditions. 

3.1 Quantization Of The Level 

In order to see the connection between the 3-dimensional Chern-Simons term and the CVE, we 
dimensionally reduce the theory on the compactified time direction [5]. To this end, we rename the 
relevant sources: 

Qi = Ai{x), (19) 
bi = Tgoi, (20) 

which transform under gauge and coordinate transformations with respective parameters a and 
as follows: 

tti^ ai + dia, (21) 
bi^bi + T (Vjco + Voe^) , (22) 

where is the diffeomorphism parameter (x'^ x^^ + e'^). As we argued in section 2, these 3- 
dimensional gauge transformations are not anomalous, even in the presence of anomalies in the 
underlying 4-dimensional theory^. The term of concern in the effective action is: 

Seff = iK J d^xeijkQidjbk = ikT J d'^xeijkaidjbk, (23) 

where in the final equality we have inserted a spurious factor of T along with a trivial time integral 
in order to make contact with the 4-dimensional 2 point function at hand. With this definition, 
the Green's function becomes: 



For coupling to gauge fields, this is only true at leading order in 1/N. 
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which leads to a contribution ot the CVE coefficient of 

a\ = kT^. (25) 

To see the quantization of the coefficient k, we put the system on a 3 torus of sides Li, L2 and 
L3. The periodicity of the manifold dictates the periodicty of the compact gauge fields a, and bi. 

Using large gauge transformations a = — and cq = we see: 

e Li T Li 

ai = ai^ (26) 

hi = hi + |-, (27) 

where e is the electric charge and the factor 2 in the time translation parameter €0 is compensating 
for the anti-periodicity of fermions in the time direction. 

We then perform a large gauge transformation on the fields with a = 27rn , with n G Z. 

The change in the action is: 

5Seff = 27mij J ezjkdjh- (28) 

X3=0 

Choosing a field configuration with non-trivial winding around the compact directions such as 

2x1 

62 = , the variation becomes: 

L1L2 

SSgff = 27rm— / £3125162 = 47rm— , (29) 
c J c 

X3=0 

which gives us the desired quantization result of k G -Z.^ 
3.2 Zero Coupling Example 

We now look at the Kaluza-Klein reduction of the action (2) at zero coupling. The fermions 
decompose on the thermal circle with periodicity /3 = ^ as: 

^{t,x) = e^^'"'^Mx), neZ + i (30) 

and we assume the background sources have no time dependence. At zero coupling the different 
Matsubara frequencies do not mix and we can analyze the partition function mode by mode. For 
the ra'th mode we have: 

S=^J(fx {ii>-n0^n + A^Jl^ - h^iTf) , (31) 



^This is half the expected result of integer n values and is solely due to the fact that the fermions are anti-periodic 
in the thermal direction as stated under equation (27). 
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with 



5,n 



Til' =^^-nr{d' - d')^n - mri^-nfi^n. (32) 



The diagram to be evaluated is given in figure la with no internal vertices. We have: 

^ J (27r)3 (p + kf + 47r2m2r2 V4r^^ ^ ^ ^ J + An^m^T^ 

where denotes contraction with 7% the spatial gamma matrices. We have: 



,,„5.* = - T / ^ „ , -"^:2'r!f. .... . (34) 



(27r)3 ((p + A;)2 + 47r2m2T2)(p2 + 47r2m2r2) ' 
where 

iV4 =i-27rmT){-7rmd^^)ps - (p + k)r^{2p + fc^p^ + (p + /c)r(-vrm(5^)(-27rmT) 

= - - ■^krPs{2p + ky + Sym in r ^ s. (35) 

Using the above and taking the limit of A; — )• 0, we have: 

,0, _^ /■ A 6,,„46'--'(27r2m2r^^fc, + ^W) 

iim6i,„C, _J y ^ (p2 + 4^2^2^2)2 • W 

This integral is of course divergent and needs to be regularized. Here we simply use the relation 
between different integrals in dimensional regularization^ . Pauli-Villars regularization gives the 
same result: 



00 

lim CiinG'''^^ =16m^T^kn f dp 
k^o J 



p2 



(p2 + 47r2m2r2)2' 


=2\m\Tkn. (37) 

And finally we write down the contribution to the CVE coefficient for each mode: 

= \m\T. (38) 

Which satisfies the quantization condition in the previous section. We evaluate the full CVE 
coefficient: 

00 00 ^ 

-\ = T Y: <m = 2T^J2m^ \-T\ (39) 

Where we have used the C- function regularization. This is the same result as [8]. In essence, this 
calculation is the same with the order of the 3-dimensional momentum integral and the Matsubara 
sum. 
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Here we just need the simplest of these relations: 



f (fl I'^r _ -A [ _dn_ 

J (277)d (Z2 + A)2 ~ d-2] (277)° 
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4 Conclusion 



We have shown in detail that, in so far as the perturbative loop expansion is valid, the zero 
chemical potential, temperature dependent part of the chiral vortical effect coefficient does not 
receive any radiative corrections from Yukawa type couplings. We have also shown that this non- 
renormalization is related to the quantized nature of the 3-dimensional effective action. 

None of this holds, however, when there are couplings to dynamical gauge fields present in the 
theory. In this case, even at large N, radiative corrections are still present and we have calculated 
the leading correction to the CVE coefficient in this limit. 

It is curious that the actual value of the CVE coefficient does not satisfy the quantization con- 
dition, despite the fact that the contribution from each Matsubara mode does. Considering the fact 
that the proof of the quantization requires invariance under both large axial gauge and large diffeo- 
morphism transformations, it seems very probable that this value would be fixed by considerations 
of global anomalies [17, 18]. This would also explain why arguments based on anomalous Ward iden- 
tities such as in [6], which are consequences of anomalies only in the infinitesimal transformations, 
would not be able to determine this coefficient. 
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